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Abstract. Radiation pressure from the absorption and scattering of starlight 
by dust grains may be a crucial feedback mechanism in starburst galaxies and the 
self-gravitating parsec-scale disks that accompany the fueling of active galactic 
nuclei. I review the case for radiation pressure in both optically-thin and highly 
optically-thick contexts. I describe the conditions for which Eddington-limited 
star formation yields a characteristic flux of ~ 10^'^ L0 kpc~^, and I discuss 
the physical circumstances for which the flux from radiation pressure supported 
disks is below or above this value. In particular, I describe the young stellar 
disk on ~ 0.1 pc scales at the Galactic Center. I argue that its bolometric 
flux at formation, ~ 10^^ L0 kpc~^, and the observed stellar mass and scale 
height imply that the disk may have been radiation pressure supported during 
formation. 



1. Introduction 

Star formation in galaxies is observed to be globally inefficient, with only a few 
percent of the available gas supply converted to stars per dynamical timescale 
(e.g., Kennicutt 1998). On smaller scales within individual galaxies star for- 
mation is similarly slow (e.g., Kennicutt et al. 2007; Bigiel et al. 2008; Leroy 
et al. 2008; Krumholz & Tan 2007). This, together with the fact that galaxies 
are marginally-stable, with Toomre's Q ~ 1 (e.g., Martin & Kennicutt 2001), 
suggests that they are self-regulated — otherwise, disks should quickly evolve 
to have Q < 1, and fragment on a single dynamical time. 

The injection of energy and momentum into the interstellar medium (ISM) 
by stellar processes ("feedback"), may be responsible for the inefficiency of star 
formation in galaxies and for their self-regulation. The feedback mechanism 
thought to stave off self-gravity is uncertain, but may be due to supernovae, stel- 
lar winds, expanding HIT regions, cosmic rays, and radiation pressure of starlight 
on dust (e.g., McKee & Ostriker 1977; Chevaher & Fransson 1984; Silk 1997; 
Wada & Norman 2001; Scoville et al. 2001; Matzner 2002; Thompson, Quataert, 
& Murray 2005 [TQM]; Li & Nakamura 2006; Socrates et al. 2007). Non-stehar 
processes can also act as a feedback mechanism by tapping the galaxy's gravi- 
tational binding energy (Sellwood & Balbus 1999). 

Both observations and theory suggest that whatever the mechanism, feed- 
back must be yet more effective on parsec scales around active galactic nuclei 
(AGN), at surface densities in excess of those probed by the densest galaxies 
that comprise the Schmidt Law. On scales larger than ~ 0.01 pc accretion disks 
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feeding bright AGN are gravitationally unstable and prone to fragmentation and 
star formation (Kolykhalov & Sunyaev 1980; Shlosman et al. 1989, 1990; Ku- 
mar 1999; Wada & Norman 2002; Goodman 2003). This problem is severe: gas 
accumulating on ~ 0.01 — fewpc scales to fuel a supermassive black hole may 
never get there, but instead fragment entirely into stars. Both rapid inflow due 
to efficient angular momentum transport (e.g., bars) and/or a reduced star for- 
mation efficiency have been proposed to help solve this "AGN fueling problem" 
(Shlosman et al. 1990; Levin 2007; TQM). 

A number of observational results also suggest an intimate connection be- 
tween AGN fueling and nuclear star formation. For example, Storchi-Bergmann 
et al. (2005) and Davies et al. (2006), (2007) find direct evidence for star for- 
mation in the central tens of parsecs of local AGN. Likewise, the young stellar 
disk on ~ 0.1 pc scales in the Galactic Center (e.g., Paumard et al. 2006; Lu et 
al. 2008) may be a vestige of accretion, providing further evidence for a close link 
between BH fueling and star formation (Levin & Beloborodov 2003; Milosavlje- 
vic &: Loeb 2004; TQM; Nayakshin et al. 2007). There is further evidence that 
the fueling of bright high-z quasars is accompanied by intense star formation 
(e.g., Lutz et al. 2007, 2008; Younger et al. 2008). 

In this contribution, I review the arguments for the importance of radiation 
pressure of starlight on dust grains as an important feedback mechanism in 
intense starburst galaxies and in the self-gravitating disks that likely attend 
AGN fueling. In I discuss some of the physics of this particular feedback 
mechanism. In ^3^ I discuss some examples of starburst systems for which the 
Eddington limit may be relevant, including local ULIRGs and the young stars 
at the Galactic Center. Section S?] provides a brief conclusion. Throughout, I 
borrow from the more extensive discussion presented in TQM. The interested 
reader is referred there for more details. 

2. Radiation Pressure Feedback 

The light from massive stars is efficiently absorbed and scattered by dust grains, 
providing a net force in the direction normal to a star-forming disk. TQM showed 
that this might be the dominant feedback mechanism in dense optically-thick 
starbursts and AGN disks (see also Scoville et al. 2001; Scoville 2003; Levin 
2007). A basic requirement for radiation pressure to be important in governing 
the dynamics is that the radiative flux approach the Eddington flux 



where Stot is the total disk surface density and (up) is the column-averaged 
flux-mean dust opacity (e.g., Elitzur &; Ivezic 2001). Note that we expect the 
gas and dust to behave as a single collisionally-coupled fluid on the scales and at 
the densities of interest (e.g., Laor & Draine 1993; Murray et al. 2005 [MQT]). 
Because the radiative flux from galaxies is proportional to the star formation 
rate per unit area (-F* = eS^c^)I3 equation ([1]) is an interesting and testable form 



~ 10 is an IMF-dependent constant; see Kennicutt (1998). 
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for the Schmidt Law: 

ec \{kf)J cm-4 V(kf)/ 

Although Stot is relatively easy to measure or estimate in normal star-forming 
galaxies and starbursts, the coupling between the radiation field and the dust, as 
expressed in equation ([T|) with (kf) is more difficult, since the medium is highly 
turbulent, inhomogeneous, and clumpy. Three simple limits help to illustrate 
the range of the effective flux- mean opacity in galactic contexts. 

2.1. Optically-Thin to UV: < 5 Mq pc"^ « 10^ g cm ^ 

When the average medium is optically-thin to the UV radiation from massive 
stars one finds that for young stellar populations (< 5Myr), standard grain size 
distributions, and Galactic dust-to-gas ratio that (kf) ~ 10'^ cm^ g"-*^. This 
limit is only applicable for galaxies with gas surface densities of < 5 Mq 
pc~^, again assuming Galactic dust-to-gas ratio. For typical numbers, 

F.,^,.0.10'L,.pc-(^^)(l2^)^ (3) 

2.2. Optically-Thick to UV, but Optically-Thin to Re-Radiated FIR: 

5M0 pc-2 < < 5000 Mq pc-2 

In normal star- forming galaxies and some starbursts, the average gas surface 
density T,g is high enough that the medium is optically-thick to the UV emission 
from massive stars, but not sufficiently large that the re-radiated FIR emission 
from dust grains is optically-thick. 

In this "single-scattering" limit (see MQT, TQM), UV photons are absorbed 
once, and then escape the system. For a homogeneous slab of gas and dust and 
with an incident UV radiation field, it is easy to show that {kf) ~ (e.g., 
TQM). Although the applicability of the homogeneous slab is highly uncertain 
in turbulent star-forming environments, the single-scattering approximation is 
useful at the order-of-magnitude level for gauging the dynamical importance of 
the radiation field. In this case. 



F,,Edd ~ 27rGStotSgC ~ 10^ L© kpc-^ ( ) { ,,,, ) . (4) 



100 Mq pC-V VlO Mq pC- 



2.3. The Optically-Thick Limit: Sg > 5000 M© pc"^ 

When the galaxy is optically-thick to the re-radiated FIR emission from dust 
grains, (kf) can be approximated by the Rosseland-mean opacity, kr{T), where 
T is the midplane temperature of the star-forming disk. Although the tem- 
perature dependence of the opacity is fairly complicated, to rough approxima- 
tion when T < 200K, ?a 2 X lO^^T^ = kqT'^, and for 200K < T < T^uh, 
kr{T) X constant ^ 5- lOcm^ g-^ (e.g.. Bell & Lin 1994; Semenov et al. 2003). 
For temperatures above the sublimation temperature of dust Tgub ~ 1500 K, the 
Rosseland mean opacity decreases markedly. 
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The Eddington Limit when T < 200 K: In this Umit, the midplane tem- 
perature is connected to the effective temperature and the radiated ffux by 
T^ ~ tTJj ~ (KoT^Sg/2)(F/(TsB)- Additionahy, if the radiation pressure bal- 
ances gravity, pr = ^ vrGSgEtot- These two equations imply that 

i^.,Edd ~ {37rcGasB/4)^^^ (Stot/Sg)'/' - lO^^ kpc-2 (^^^^/^^f/^ , (5) 

where kq = 2 x 10~^ cm^ g~^ has been assumed. For typical values of e 
(eq. [2|), equation ^ predicts a star formation rate surface density of approx- 
imately ~ 10^ Mq yr~^ kpc~^ required to support the disk with radiation 
pressure on dust. Note that because we expect starbursts that meet the criteria 
used to derive equation ([5|) to be largely gas-dominated, Stot is a reason- 

able first approximation. Hence, starburst galaxies with T,g > 5000 pc^^ 
should attain a characteristic flux of ~ 10^^ Lq kpc~^. 

The primary assumptions made in deriving equation ([5]) are that the medium 
is optically-thick and that the midplane temperature is less than 200 K. For 
example, a massive starburst with a gas reservoir of - 10^° M0 distributed on 
scales < 1 kpc should satisfy both criteria over approximately 1 decade in radius, 
down to scales of < 100 pc, depending on the radial distribution of the stellar 
and gas mass. If the surface density increases at smaller radii, the midplane 
temperature will as well and the assumption that T < 200 K will be violated. 
Similarly, on somewhat larger scales, the medium is optically-thin to the re- 
radiated FIR emission and equation becomes applicable. 

The Eddington Limit when T > 200 K: If the equilibrium radiative flux 
required for hydrostatic equilibrium implies 200 K < T < T^^h^ then FEdd 
exceeds the value given in equation ([5]). In this case, we can approximate 
kr{T > 200 K) ~ 5 - lOcm^ g"! and 

F.,Edd~10 Lekpc j Uo^MePC-^ j' 

where I have scaled the total surface density for parameters characteristic of 
~ 0.1 — l.Opc self-gravitating disks around AGN (see ^3. 2. p . 

2.4. Stability 

Equation ([5]) implies that if an optically-thick starburst with T < 200 K reaches 
a characteristic flux of ~ 10^^ Lq kpc~^, then it can be maintained in hydro- 
static equilibrium in the sense that radiation pressure balances the self-gravity 
of the disk. However, because the diffusion timescale for radiation is very short 
compared to the dynamical timescale, the disk is unstable. In other words, ra- 
diation pressure cannot stave off the local Jeans instability on the scale of the 
disk gas scale- height because diffusion is rapid (see Thompson 2008). This ef- 
fect should keep the disk highly turbulent, with the magnitude of the resulting 
turbulent energy density {p5v^) set by the magnitude of the radiation energy 
density. In this way, hydrostatic equilibrium would be maintained statistically. 
As discussed in Thompson (2008), turbulence may also generate the ~ fewmG 
magnetic field strengths needed in the densest starbursts to explain the fact that 
these systems lie on the FIR-radio correlation (Thompson et al. 2006). 




Figure 1. Left: Flux as a function of radius, as computed by dividing the 
FIR luminosity by the radio continuum area, for the local ULIRGs in Condon 
et al. (1991; filled squares) compared with a set of theoretical models from 
TQM using equations ^ and ^ (solid & dashed lines). The solid lines 
have constant gas fraction of 0.1, 0.3, and 1.0 from lowest to highest. The 
dashed line shows a model with a deeper stellar potential. Right: Histogram 
showing number of ULIRGs as a function of flux. The peak at ~ 10^'^ L0 
kpc~^ is consistent with the prediction of Eddington-limited star formation 
in radiation pressure supported self-gravitating disks. Plots from TQM. 



3. Some Examples 

3.1. Local & High-z Starbursts 

The flux predicted by equations ([5]) and ([6]) can be compared with observations 
of local and high-z ultra-luminous infrared galaxies. Figure [T] {left) shows a 
comparison between the observed fluxes of local ULIRGs and a set of models 
with increasing constant gas fraction from TQM. The right panel displays the 
same data as a histogram, and indicates that ULIRGs may in fact exhibit a 
characteristic flux of ~ 10^^ Lq kpc~^. It is also interesting in this context 
to note that Forster-Schreiber et al. (2003) find that the ~ lOMyr starburst 
required to explain the observed stellar population in the central starburst of 
M82 also attained a star formation rate surface density of ~ 10^ Mq yr~^ 
kpc~^ (see their Fig. 13). 

Finally, observations of submillimeter galaxies at high-z are also in fair 
agreement with the predictions of TQM; Younger et al. (2008), Walter et al. (2009), 
and Riechers et al. (2009) argue that the luminosities of the systems they observe 
indicate that radiation pressure on dust may be dynamically important. 

3.2. The Galactic Center & The Fueling of AGN 

A number of works have shown that there is disk of young stars around Sag A* at 
the Galactic Center with total stellar mass ~ 2 X 10^ Mq, radius r ~ 0.1 pc, 
scale height of /i ~ O.lr, and an average age of ~ 6Myr (Levin & Beloborodov 
2003; Genzel et al. 2003; Paumard et al. 2006; Lu et al. 2009). Although the 
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physical origin of these stars is still unclear, it is possible that the stars formed 
in situ in a gaseous disk (Levin & Beloborodov 2003). This star formation 
episode may have been accompanied by accretion onto the central black hole, 
perhaps producing an AGN epoch in the Galaxy (e.g., TQM; Nayakshin &: 
Cuadra 2005; Levin 2007). Importantly, its structure may have been determined 
by the radiation pressure of starlight on dust grains. 

A simple estimate for the minimum star formation rate needed to produce 
the stars that we see now is M™™ ~ Mass/Age ~ 3 x 10~^ Mq yi~^, cor- 
responding to a star formation rate surface density of S™™ ~ M™™/7rr^ ^ 
10^ M0 yr-i kpc-2, roughly 100 times the surface density of star formation for 
the central 100 pc of Arp 220's western nucleus. Note that the ratio L^/M^ for 
a few-Myr old stellar population is 1^1500 ~ 1500 Lq/Mq, which allows for an 
order-of-magnitude estimate for the luminosity and flux of the stars at forma- 
tion: ~ 3 X lO'^^isooL© and ~ 10^^^i5oo kpc~^, respectively. 

It is interesting to consider the physical conditions of the disk in forma- 
tion using the observed stellar mass density. We can write down an approxi- 
mate lower limit to the total gas mass required to build the stars we see today 
by equating the gas mass and the stellar mass: Mg ~ M^. This implies a 
minimum volumetric and surface gas density of n ~ Mg/{2TTr^h) ~ 10^ cm 
and ~ 2ph ~ 120 g cm~^ '--^ 6 x 10^ Mq pc~^, respectively, assuming that 
the stars formed with h ~ O.lr. Note that this estimate for the gas density 
of the disk during formation implies an average dynamical timescale of order 
tdyn ~ {Gp)^^/'^ ~ 3000 yr, which implies that the maximum star formation rate 
exceeds M™" by a factor of ~ 2000: Aff^^ ~ GMoyr'^ 

Given these approximate parameters for the disk in formation, we can ask if 
it is physically possible for a given mechanism to support the gas in vertical hy- 
drostatic equilibrium with scale height of order h. The effective surface tempera- 
ture of the disk is T^q ^ (F/o'sb)^''^ ~ 260 K and the total average optical depth 
in the vertical direction is r ~ 600^10^9/10. ir where kiq = k_r(T > 200 K)/10 cm^ 
g~-^. Note that this implies a midplane temperature of ~ 1300 K, close to — but 
somewhat below — the sublimation temperature of dust. Given T and n, the 
ratio of the radiation pressure to the gas pressure is of order Pr/Pg ~ 20. 

Hydrostatic equilibrium in the vertical direction requires that the radiated 
flux equal the Eddington flux. Because the midplane temperature exceeds ~ 
200 K, the estimate of equation ^ is appropriate and the ratio of the Eddington 
flux to the flux from starlight during formation is 

FEdd 27rGc -it-i 

—E' — ^ ~ 1 ^10?1500' (') 

where ~ lO^^^isoo L© kpc^^. The fact that equation ([7]) yields a number close 
to unity demonstrates the importance of radiation pressure in this environment 
for maintaining the disk in vertical hydrostatic equilibrium. Although there 
is evidence for a top-heavy IMF (Nayakshin & Sunyaev 2005; Nayakshin et 
al. 2006), there is relatively little uncertainty in ^ in the above expression relative 
to the assumption that the disk of observed stars can be treated as a disk 
of gas with approximately the same geometry. It should be noted again that 
the uncertainty in k/j is dominated by the assumption of Galactic dust-to-gas 
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Figure 2. Predicted star formation rate (Af^; solid lines) and gas accretion 
rate (M; dashed lines) as a function of radius for conditions appropriate to 
the Galactic Center, using the model of Eddington-limited star formation in 
TQM and discussed in ij2] The different line sets indicate different input gas 
accretion rate at a fiducial outer radius of 3 pc, varying from Mout = 0.015 
yr~^ to 0.015 Mq yr~^. Solutions with higher Mout produce a peak in S^, of 
~ 10^ - 10*^ M0 yr-i kpc-2 on 0.1 pc scales. 



ratio and Solar metallicity. If anything, k\q may be an underestimate of the 
true opacity. Nevertheless, equation ([7]) highlights the importance of radiation 
pressure in extreme star formation environments and may indicate that the 
formation of the stars in the Galactic center region was in fact Eddington-limited 
(see TQM & Levin 2007 for more details). 

Equation ([7|) also shows that radiation pressure supported disks may vastly 
exceed a bolometric flux of ~ 10^^ Lq kpc~^ if the midplane temperature of the 
disk is larger than 200 K (eq. [6l). In the case of the Galactic Center, the flux 
needed for radiation pressure support is of order ~ 10^^ kpc~^, and this is 
in good agreement with the inferred flux from the stars at the Galactic Center. 
This is strong evidence that radiation pressure may be paramount in the self- 
gravitating disks that accompany AGN fueling. Indeed, in a detailed study of 
star formation in local Seyfert nuclei, Davies et al. (2007) find evidence that star 
formation reaches the characteristic flux of equation ([5]) on ~ 1 — lOpc scales. 
More evidence (besides the Galactic Center) that the flux from young stars can 
exceed ~ 10^^ Lq kpc~^ on very small scales would be useful in constraining 
detailed models of AGN fueling (e.g.. Levin 2007; Nayakshin et al. 2007). 

As an example of a suite of such models. Figure [2] shows the average star 
formation rate as a function of radius from TQM (solid lines) for parameters ap- 
propriate to the Galactic Center, with a central black hole mass of ~ 4 x 10^ Mq 
and a bulge velocity dispersion of ~ 75 km s~^. For a given input accretion rate 
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at 3pc (Mout), the disk forms stars at the rate required by the Eddington hmit. 
The dashed hnes show the gas accretion rate, which declines as a function of 
radius towards the black hole as a result of star formation. Models that fuel 
the black hole with ~ 10~^ Mq yr"^ reach ~ 10^ - 10*^ M© yr~^ kpc~^ on 
0.1 pc scales, producing a sharp peak in the resulting stellar surface density after 
accretion abates (see Appendix A in TQM). 

Other Sources of Feedback: In addition to radiation pressure, it is also worth 
considering other sources of pressure support thought to be important in star- 
forming environments on galactic scales: gas, cosmic rays, and magnetic fields 
(e.g., Boulares & Cox 1990). The gas temperature required to support the disk 
is Tg ~ irGT^l/nkB ~ 2 X 10^ K. However, if the medium is optically-thick, 

then Pr/Pg ~ lO^T^ng ^, where T4 = T/IO^K. Additionally, cooling timescale 
arguments and the total energy required from the stellar or AGN radiation 
field show that Tg is difficult to maintain (although, sec the discussion of pg 
in Appendix A of TQM). The magnetic field strength required for hydrostatic 
equilibrium is given by B^/Stt ~ ttGS^, which implies B ~ 0.3 G for fiducial 
parameters. For cosmic rays, I assume that a fraction / of the total holometric 
output from the stars goes into primary cosmic ray protons: £^CR ~ Z-^*; as 
perhaps might be provided by stellar winds. The cooling timescale is dominated 
by inelastic proton-proton scattering, with a typical cooling timescale for few- 
GeV protons of ~ O.OSn-g^yr. This gives a maximum upper limit to the 
equilibrium midplane CR pressure of PcR ~ /i'*i7r/ (3F) ~ 3 x 10~®/ ergs cm~^, 
which is very small with respect to the total pressure required for hydrostatic 
equilibrium, ~ 3 x 10~^ergs cm~^. I conclude that gas pressure and cosmic ray 
pressure are both likely sub-dominant with respect to radiation pressure in this 
context, and that magnetic fields of strength approaching ~ G are required to 
contribute significantly to the pressure budget. 

Supernovae are similarly unlikely to have supported the disk, given the 
youth of the observed stellar population. 



4. Summary &; Conclusions 

Radiation pressure associated with the absorption and scattering of starlight by 
dust grains in rapidly star-forming environments is likely to be an important 
feedback process, and may be responsible for the regulation of star formation in 
some circumstances on galactic scales. Here, I have reiterated and clarified ar- 
guments in the starburst and AGN contexts developed in TQM. However, much 
more work needs to be done to assess this mechanism in full. In particular, it 
is crucial to understand the coupling between the radiation field produced by 
spatially clustered massive stars and the highly inhomogeneous and turbulent 
ISM of starbursts (e.g., Krumholz & Matzner 2009; Murray et al. 2009). This 
remains a primary avenue of future investigation. 
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